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ABSTRACT 

We characterize a cosmic rest frame in which the monopole variation of the spherically 
averaged nonlinear Hubble expansion is most uniform, under arbitrary local Lorentz 
boosts of the central observer. Using the COMPOSITE sample of 4534 galaxies, we 
identify a degenerate set of candidate minimum nonlinear variation frames, which 
includes the rest frame of the Local Group (LG) of galaxies, but excludes the standard 
Cosmic Microwave Background (CMB) frame. Candidate rest frames defined by a 
boost from the LG frame close to the plane of the galaxy have a statistical likelihood 
similar to the LG frame. This may result from a lack of constraining data in the Zone 
of Avoidance. We extend our analysis to the Cosmicflows-2 (CF2) sample of 8162 
galaxies. While the signature of a systematic boost offset between the GMB and LG 
frame averages is still detected, the spherically averaged nonlinear expansion variation 
in all rest frames is significantly larger in the GF2 sample than would be reasonably 
expected. We trace this to the GF2 distances being reported without a correction for 
inhomogeneous distribution Malmquist bias. Systematic differences in the inclusion of 
the large SFI++ subsample into the COMPOSITE and CF2 catalogues are analysed. 
Our results highlight the importance of a careful treatment of Malmquist biases for 
future peculiar velocities studies, including tests of the hypothesis of Wiltshire et al. 
(2013) [Phys. Rev. D, 88, 083529] that a significant fraction of the CMB temperature 
dipole may be nonkinematic in origin. 
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1 INTRODUCTION 

Although the Universe is spatially homogeneous in some sta¬ 
tistical sense, at the present epoch it exhibits a complex hi¬ 
erarchical structure, with galaxy clusters forming knots, fila¬ 
ments and sheets that thread and surround voids, in a com¬ 
plex cosmic web (Forero-Romero et al. 2009; Bilicki et al. 
2014; Einasto 2014). Deviations from homogeneity are con¬ 
ventionally treated in the framework of peculiar velocities, 
by which the mean redshift, z, and luminosity distance, r, 
of a galaxy cluster are converted to a peculiar velocity at 
small redshifts according to 

Upec = cz - H^r (1) 

where c is the speed of light and = 

100 h km sec~^ Mpc“^ is the Hubble constant. 

While large galaxy surveys reveal that statisti¬ 
cal homogeneity emerges on scales 70 - 100/i~^Mpc 

* E-mail: j.mckayl4@imperial.ac.uk 
f E-mail: david.wiltshire@canterbury.ac.nz 


(Hogg et al. 2005; Scrimgeour et al. 2012) as determined 
by the two-point galaxy correlation function, in the 
standard peculiar velocities framework relatively large 
bulk flows are seen on scales larger than this rel¬ 
ative to our own location. The amplitude of such 
flows, and their consistency with the standard Lambda 
Cold Dark Matter (ACDM) cosmology, are a mat¬ 
ter of ongoing debate (Watkins, Feldman & Hudson 2009; 
Kashlinsky et al. 2010; Nusser & Davis 2011; Turnbull et al. 
2012; Lavaux et al. 2012; Ma & Scott 2013; Atrio-Barandela 
2013; Ade et al. 2014a; Hoffman, Courtois & Tully 2015; 
Garrick et al. 2015). While it may be possible to achieve 
convergence of such bulk flows to the CMB frame within 
current uncertainties based on cosmic variance and A-body 
simulations in the ACDM model, given the deep mystery as¬ 
sociated with the nature of dark energy, it is important that 
alternative cosmological models are investigated. Such alter¬ 
natives (Buchert 2000, 2008; Wiltshire 2007a,b, 2009) may 
also demand a treatment of cosmic expansion which includes 
effects of inhomogeneities that fall outside the peculiar ve- 
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locity framework, with testable observational consequences 
(Li & Schwarz 2008; Wiltshire et al. 2013). 

The peculiar velocity framework makes a strong ge¬ 
ometrical assumption over and above what is demanded 
by general relativity. In particular, the quantity Upec de¬ 
fined by (1) only has the physical characteristics of a ve¬ 
locity if one implicitly assumes the spatial geometry on all 
scales larger than those of bound systems is exactly de¬ 
scribed by a homogeneous isotropic Friedmann-Lemaitre- 
Robertson-Walker (FLRW) model with a single cosmic scale 
factor, a{t), whose derivative dehnes a single global Hubble 
constant, Hp = a/a\^^. Deviations from the uniform ex¬ 
pansion are then ascribed to Lorentz boosts of each galaxy 
cluster with respect to the spatial hypersurfaces of average 
homogeneity, at each cluster location. 

It is a consequence of general relativity, however, that 
inhomogeneous matter distributions generally give rise to a 
differential expansion of space that cannot be reduced to a 
single uniform expansion plus local boosts. This is a feature 
of general exact solutions to the cosmological Einstein equa¬ 
tions, such as the Lemaitre-Tolman-Bondi (LTB) (Lemaitre 
1933; Tolman 1934; Bondi 1947) and Szekeres (1975) models. 
Any definition of the expansion rate in such models depends 
on the spatial scale relative to that of the inhomogeneities. 
Although one can define scale dependent Hubble parameters 
for specihc exact solutions - for example, given the spherical 
symmetry of the LTB model - the actual cosmic web is suf- 
hciently complex that in reality one must deal with spatial 
or null cone averages in general relativity. 

In recent work Wiltshire et al. (2013) examined the 
variation^ of the Hubble expansion from a fresh perspec¬ 
tive, by generalizing the earlier approaches of Li & Schwarz 
(2008) and McClure & Dyer (2007). In particular, given 
that there is a notion of statistical homogeneity on large 
(3> 100 h“^Mpc) scales, then an average expansion law char¬ 
acterized by a single asymptotic Hubble constant, H^, is ap¬ 
plicable on such scales. However, from the first principles of 
general relativity one should make no geometrical assump¬ 
tions about cosmic expansion below the statistical homo¬ 
geneity scale, where cosmic expansion is nonlinear. One can 
nonetheless perform radial and angular averages of the dis¬ 
tance versus redshift of a large sample of galaxies in spherical 
shells, and compare the results with the asymptotic Hubble 
constant in order to quantify the nonlinear variation of the 
Hubble expansion. 

Wiltshire et al. (2013) conducted such an analysis 
on the COMPOSITE sample of 4534 cluster, group 
and galaxy distances (Watkins, Eeldman & Hudson 2009; 
Feldman, Watkins & Hudson 2010), with the following re¬ 
sults: 

• A linear Hubble law with a spherically averaged Hubble 
constant, Ha, which is statistically indistinguishable from 
the asymptotic Hubble constant, H^, is found to emerge in 
independent radial shells with mean distances in the range 
fa > 70 h“^Mpc. 

^ Wiltshire et al. (2013) used the terminology “Hubble flow vari¬ 
ance”, where variance was understood in the same loose sense 
as in the terminology “cosmic variance”. Here we use the more 
generic word “variation”, to avoid any confusion with the square 
of a standard deviation. 


• On scales r < 65/i“^Mpc the spherically averaged value, 
Ha, in independent shells is greater than the asymptotic 
value, Hq. However, the difference is significantly larger in 
the standard rest frame of the CMB radiation when com¬ 
pared to either the rest frame of the Local Group (LG) of 
galaxies or the Local Sheet (LS). In other words, the spheri¬ 
cally averaged Hubble expansion is more uniform in the LG 
rest frame than in the GMB rest frame, with very strong 
Bayesian evidence, InB ^ 5. The uniformity of expansion 
in the LG and LS frames is statistically indistinguishable. 

• By a variety of angular tests, the residual variation of the 
spherical (monopole) Hubble expansion is found to be cor¬ 
related with structures in the range 32 - 62 h“^Mpc, which 
give a Hubble expansion dipole with a markedly different 
character in the GMB and LG frames. 

• A skymap of angular variation of the Hubble expansion in 
the LG frame, constructed by Gaussian window averaging 
(McGlure & Dyer 2007), has a very strong dipole. The angu¬ 
lar expansion skymap has a correlation coefficient of —0.92 
with the residual CMB temperature dipole in the LG frame. 

• On scales > 80h“^Mpc the magnitude of the dipole vari¬ 
ation in the Hubble expansion is less in the LG frame than 
in the CMB frame. While the LG frame dipole is statisti¬ 
cally consistent with zero in most individual outer shells, 
by contrast the CMB frame dipole magnitude rises to a 
residual level. Consequently the large bulk flow reported by 
Watkins et al. (2009) may be partly due to a choice of rest 
frame. 

The Erst of the results above is consistent with other ob¬ 
servations which find that a notion of statistical homogene¬ 
ity emerges at scales of order 70 - 100 h“^Mpc (Hogg et al. 
2005; Scrimgeour et al. 2012). Furthermore, the fact that 
Ha > Hq on the < 65/i“^Mpc scales on which the Hubble 
expansion is nonlinear agrees well with the observation that 
the largest typical^ structures in the late epoch Universe are 
voids of diameter 30h“^Mpc (Hoyle & Vogeley 2002, 2004). 
Pan et al. (2011) found that voids occupy 62% of the vol¬ 
ume studied in the 7th release of the Sloan Digital Sky Sur¬ 
vey (Abazajian et al. 2009). If one constructs averages in 
spherical shells, as in Fig. 1, then once the shells are a few 
times larger than the diameter of the largest typical non¬ 
linear structures, a well-defined asymptotic average, H^, is 
obtained which does not change when shells are further en¬ 
larged. When shells are 1-2 times the diameter of the 
typical nonlinear voids, however, a variation in expansion 
rate is seen and since the faster expanding voids dominate 
by volume then the average, Ha, is increased relative to H^. 

While the results just discussed may not be surprising 
to those familiar with the statistics of the cosmic web, the 
remaining results of Wiltshire et al. (2013) listed above are 
not at all expected in the conventional approach to pecu¬ 
liar velocities in observational cosmology. In particular, the 
cosmic rest frame should be the one in which the variation 
of cosmic expansion is a minimum in some sense. The re¬ 
sults of Wiltshire et al. (2013) show that such a frame of 

^ Larger structures such as the 320/i“^Mpc long Sloan Great 
Wall (Gott et al. 2005) and the 350/r“^Mpc long Large Quasar 
Group (Clowes et al. 2013) are known, but it is arguable that 
these are not typical. 
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Figure 1. Schematic diagram of spherical averaging (from 
Wiltshire et al. (2013)). The universe is an ensemble of filaments, 
walls and voids, the largest typical nonlinear structures being 
voids of diameter ~30/i“^Mpc (Hoyle & Vogeley 2002, 2004; 
Pan et al. 2011). If one averages cz/r in spherical shells (dot¬ 
ted lines) about any point then once the shells are a few times 
larger than the typical nonlinear structures, an average Hubble 
law with small statistical scatter is obtained (e.g., for the longest 
null geodesics depicted by lines converging to the central ob¬ 
server). However, there are considerable deviations for shells on 
scales comparable to the typical nonlinear structures. This ap¬ 
proach makes no assumptions about spatial geometry of the uni¬ 
verse on the scales in question, and thus is model independent. 


minimum nonlinear Hubble expansion variation is not the 
conventional CMB frame. 

One important consequence is that a substantial frac¬ 
tion of the CMB dipole usually attributed to a local boost 
of the Local Group at 635 ±38 km sec“^ (Tully et al. 2008) 
in the direction {£,b) = (276.4°, 29.3°) would be due to a 
nonkinematic differential expansion of space. It has been 
independently shown that the radio galaxy number count 
dipole appears to be nonkinematic, at the 99.5% confidence 
level, using the NRAO VLA Sky Survey (Rubart & Schwarz 
2013). 

Wiltshire et al. (2013) estimate that a 0.5% anisotropy 
in the distance-redshift relation on <65/i~^Mpc scales 
would be required to achieve the observed properties of the 
CMB dipole. Furthermore, using ray-tracing simulations on 
an exact Szekeres solution of Einstein’s equations which ex¬ 
hibits differential expansion on such scales, it is found that 
it is possible to produce a Hubble expansion dipole which is 
very close to that of the COMPOSITE sample, while simul¬ 
taneously generating CMB anisotropies that both account 
for the residual temperature dipole in the LG frame and 
are consistent with the amplitude of higher order multipoles 
(Bolejko, Nazer & Wiltshire 2015). 

This may potentially provide an explanation for the an¬ 
gular scale dependence observed in attempts to measure the 
boost to the CMB frame using the effects of special rela¬ 
tivistic aberration and modulation on the CMB anisotropy 
spectrum (Aghanim et al. 2014). Since subtraction of a kine¬ 
matic dipole is a key step in the map making process, this 
may also have a direct impact on various large angle anoma¬ 
lies in the CMB anisotropy spectrum (Ade et al. 2014b). 
Furthermore, given that a transformation to the CMB frame 
is a step that is taken in many observational procedures 
without question, the possibility that a substantial fraction 


of the CMB dipole may be nonkinematic could have a large 
impact in other areas of observational cosmology. 

Given the potential importance of such a result, it is 
important to try to characterize the frame of minimum non¬ 
linear Hubble expansion variation in purely observational 
terms. Wiltshire et al. (2013) compared the LG and LS 
frames with that of the GMB, motivated by the starting 
point that a frame close to the LG frame would be the nat¬ 
ural standard of rest according to the “Gosmological Equiv¬ 
alence Principle” (Wiltshire 2008) which underlies the ap¬ 
proach of Wiltshire (2007a,b, 2009) to the averaging prob¬ 
lem in inhomogeneous cosmology (Ellis 1984; Buchert 2000, 
2008; Wiltshire 2011; Nazer & Wiltshire 2015). 

This paper asks the question - if we make no assump¬ 
tions about our own standard of rest, but make arbitrary 
boosts at our own location - can we determine in a model- 
independent fashion whether there is a rest frame in which 
the nonlinear variation of the Hubble expansion is mini¬ 
mized? We will consider all scales over which cosmic expan¬ 
sion is nonlinear and transitions to linearity, and the ques¬ 
tion of how such a transition impacts the characterization 
of a minimum variation frame. 

While our principal results are determined from an 
analysis of the GOMPOSITE sample (Watkins et al. 2009; 
Feldman et al. 2010), we have also considered the re¬ 
cently published Cosmicflows-2 (CF2) sample of Tully et al. 
(2013). We find that the issue of the different treatments of 
Malmquist bias in the two datasets at present prevents as 
detailed analysis as we perform for the GOMPOSITE sam¬ 
ple. However, our investigation highlights how the implicit 
assumption of a FLRW expansion law below the scale of sta¬ 
tistical homogeneity via (1) subtly influences the manner in 
which such biases are treated in practice, and raises concerns 
about remaining distribution biases in the CF2 distances. 


2 METHODOLOGY 

In this paper we will use the methodology of Wiltshire et al. 
(2013), who considered both radial and angular aver¬ 
ages of the Hubble expansion in the GOMPOSITE sam¬ 
ple, adapting techniques that Li & Schwarz (2008) and 
McGIure & Dyer (2007) had previously applied to radial and 
angular averages respectively in the much smaller Hubble 
Telescope Key Project dataset (Freedman et al. 2001). We 
are primarily interested in the radial spherical averages, as 
it was by this method that Wiltshire et al. (2013) found the 
most decisive evidence that the Hubble expansion is signif¬ 
icantly more uniform in the LG or LS frames, as compared 
to the standard GMB rest frame. 

We begin by repeating the analysis of Wiltshire et al. 
(2013), but without assuming that our own standard of rest 
is known. In particular, we will make arbitrary local boosts 
of the central observer, and perform numerical studies on 
the COMPOSITE sample to determine what cosmic rest 
frame(s) are singled out by minimizing the variation in the 
spherically averaged nonlinear Hubble expansion. 

An arbitrary boost of the central observer has the effect 
of transforming the inferred redshift of individual sources we 
observe from values, Zi, in any one rest frame to new values 
given by 

cz'i = CZi + V cos 4>i (2) 
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where v is the boost magnitude at our origin, and <j)i is the 
angle between each data point and the boost direction. We 
stress that in our nonlinear framework the boost velocity v of 
the central observer is a distinct quantity from the peculiar 
velocities of individual sources^, Upec, in the linear peculiar 
velocities framework (1). 

Following Wiltshire et al. (2013), we determine the best 
fit linear Hubble law - even in the regime in which the av¬ 
erage expansion is nonlinear - by standard linear regression 
in independent radial shells, minimizing the quantity 

Ns 

X^s ='^[K^{ri - cZi/H)]'^ , (3) 

i=l 


with respect to H, where Zi, Vi and at are respectively the 
redshift, distance and distance uncertainty of each object, 
and Ns is the number of data points in a given shell, s. 
This method, with distance as a function of redshift, is cho¬ 
sen because all uncertainties in the COMPOSITE sample 
have been included as distance uncertainties. Any correc¬ 
tions that would be required due to noise arising from pe¬ 
culiar motion within galaxy clusters have been accounted 
for by assigning a distance and uncertainty to the cluster 
itself rather than individual galaxies (Watkins et al. 2009). 
The values of Zi are taken to be exact, whereas the dis¬ 
tance uncertainties are large, approximately 15% for most 
COMPOSITE sample objects. Fortunately, the large num¬ 
ber of data points within each shell nonetheless lead to sta¬ 
tistically significant results. 

The value of the Hubble constant in the sth shell is then 
determined as 


Hs 




( 4 ) 


while its weighted average luminosity distance is 


( 5 ) 


The radial averages are computed in two different shell 
configurations of 11 shells, for luminosity distances with 
Ts < r ^ rs+i. Both configurations use shells of width 
12.5/i“^Mpc in most cases, but start from a different in¬ 
nermost shell cutoff of either 2h~^Mpc or 6.25h“^Mpc. 
The two shell configurations are labelled using unprimed 
and primed integers respectively, as given in Table 1 of 
Wiltshire et al. (2013). Given less data at large distances, 
the shells 10 and 10' are taken to be wider than the rest 
so as to include approximately the same number of data 
points as the inner shells. They both have an outer cutoff at 
156.25 /i“^Mpc. The outermost shell 11 has an outer bound 
given by the largest distance in the sample, and is the same 
in both configurations, to provide an anchor for the asymp¬ 
totic value of the Hubble constant, H^. 

The uncertainty in (4) is taken as 


a 


S 



+ d? 


( 6 ) 


^ In particular, our minimum nonlinear Hubble expansion vari¬ 
ation frame should not be confused with the terminology 
of Watkins et al. (2009); Feldman et al. (2010); Agarwal et al. 
(2012), who use the term “minimum variance” with respect to 
weights of peculiar velocities in bulk flow studies. 


where 




and 


E 


iczW 


3/2 


E 


CZiTi 


( 7 ) 

( 8 ) 


A zero point uncertainty (7) is added in quadrature to the 
standard uncertainty (8) determined through error propa¬ 
gation in (4), to yield the total uncertainty (6) for the Hub¬ 
ble constant in each shell. The zero point uncertainty arises 
from the fact that a linear Hubble law must necessarily pass 
through the origin. However, there is an uncertainty in de¬ 
termining the origin, due to a 20 km sec“^ uncertainty in 
the heliocentric peculiar velocity of the LG (and LS) frames 
(Tully et al. 2008) and a 0.4% uncertainty in the magnitude 
of the CMB dipole (Fixsen et al. 1996), which combine to 
give (To = 0.201 h“^Mpc (Wiltshire et al. 2013). For each 
shell, the weighted zero point uncertainty (7) is the uncer¬ 
tainty in the linear Hubble law for a shell with mean dis¬ 
tance fs due to the uncertainty in the origin alone. This 
uncertainty is more significant for shells with a small mean 
distance compared to those at a larger distance, due to the 
shorter lever arm of the linear Hubble law for shells with 
smaller mean distances. 

As a measure of the difference in the local Hubble 
expansion from its asymptotic value, we use the quantity 
(Li & Schwarz 2008; Wiltshire et al. 2013) 


■ iW-H,) 
Ho 


( 9 ) 


where is the mean asymptotic value of the Hubble 
constant. In our case, Lfg, and its uncertainty are calcu¬ 
lated from the data points at distances r > 156.25 h“^Mpc 
(shell 11 in both unprimed and primed configurations). 
Wiltshire et al. (2013) choose the distance scale for the out¬ 
ermost shell as one which is statistically reliable for the 
COMPOSITE data, while including only data at distances 
larger than the baryon acoustic oscillation scale, the largest 
scale on which one might expect to see the effects of in¬ 
homogeneity on the local Hubble expansion. They verify 
that a linear Hubble law with a very high goodness of fit is 
found for r > 156.25 /i“^Mpc, with a Hubble constant = 
(100.1 ± 1.7)/ikm sec“^ Mpc“^ in the CMB frame or L/g = 
(101.0 ± 1.7)hkm sec“^ Mpc“^ in the LG/LS frames, con¬ 
sistent with the normalization Ho ~ 100 km sec“^ Mpc“^ 
used in the COMPOSITE sample. 


3 MONOPOLE EXPANSION VARIATION DUE 
TO SYSTEMATIC BOOST OFFSETS 

Wiltshire et al. (2013) propose that the much larger 
monopole variation of the linear Hubble law observed in the 
CMB frame, as compared to the LG/LS frames, may have 
a systematic origin. This arises from the nonlinear depen¬ 
dence of Ha in (4) on the individual czi values, which change 
when performing a local Lorentz boost to a different frame. 
By calculating the result of an arbitrary boost on the indi¬ 
vidual Ha values Wiltshire et al. (2013) obtain an explicit 
form for this systematic variation. 
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Consider redshifts, Zi, observed in a frame of reference 
in which the variation of the spherically averaged Hubble 
expansion is minimized. An arbitrary local boost (2) of the 
central observer leads to inferred redshifts, z'i, in the new 
frame given by 


law to the observed data. As there is a correlated uncertainty 
in both the independent and dependent variables, and 

AHa respectively, a standard least squares method is not ap¬ 
propriate. Instead we use a total least squares ht, or “error 
in variables” method with a model of the form 


czi cz'i = c (7 — 1) -I- 7 [czi -b V cos 0i(l -I- Zi)] 

CZi + V cos + Zi) + O (10) 

CS CZi + V cos (l>i + O , (11) 


where 7 = [l — w^/c^] We will adopt the widely used 
Newtonian velocity addition approximation (11), which is 
also assumed in (1). This ignores the 0{v/c) correction, 
which is at most 0.5% for the boosts considered here, at 
least one order of magnitude smaller than typical distance 
uncertainties. This results in the changes (czi)^ (cz))^ — 
(czi)^ -b 2cZiV cos (j)i +D cos^ (j>i in the numerator of (4), and 
cZiVi cZiVi -b ViV cos (f>i in the denominator. The linear 
contributions to the transformed quantities in the denomi¬ 
nator and numerator of (4) should be approximately self¬ 
cancelling for spherical averages of data which is uniformly 
distributed on the celestial sphere in each shell: on average 
each positive contribution from the term v cos (fn will cancel 
with a negative contribution from a data point on the oppo¬ 
site side of the sky"*. With such a cancellation assumed we 
are left with the difference 


AHa =H'a-Ha 


where 





({v cos 

lit ) 

2^0 i^i)s 




( 12 ) 


(13) 


is a weighted average of the squared luminosity distance in 
each shell. The second line of (12) follows by assuming that: 
(i) ({v cos 4>i)‘^) = is roughly constant from 

shell to shell; (ii) the leading order linear Hubble approxi¬ 
mation czi ~ Hf^ri is made in the denominator. The uncer¬ 
tainty in (13) is given by 


(T 





(14) 


Our first goal is to verify that the difference in the spher¬ 
ically averaged Hubble expansion between the LG and CMB 
frames of reference is statistically consistent with a system¬ 
atic variation of the form (12). To achieve this we fit a power 


^ The lack of data in the Zone of Avoidance does not pose a 
problem as this absence of data is symmetrical on opposite sides 
of the sky. A rough self-cancellation of the linear contributions 
would only be invalid when one side of the sky has a significant 
lack of data as compared to the opposite side of the sky. The 
COMPOSITE sample does indeed have sufficient sky coverage to 
satisfy this requirement, with the exception of the first unprimed 
shell 1, with 2 < r ^ 12.5/i“^Mpc (Wiltshire et al. 2013), which 
is excluded in the data analysis. 


AHa = A{{r^)f. (15) 

In comparison with (12), we expect p « —1 and A « 
/{2Hq). The details of this method are presented in the 
Appendix. 

Carrying out this analysis, we do indeed find a difference 
consistent with a systematic boost offset between the LG 
and CMB frames of reference. Systematic uncertainties arise 
in the choice of shell boundaries. Considering only primed 
shells gives a value of p = —1.01 ± 0.27. If we take the un¬ 
primed shells then we obtain p = —0.79 ± 0.16 if shell 1, 
with 2 < r ^ 12.5 h“^Mpc, is included and p = — 0.77T0.35 
if this first shell - which may have insufficient sky coverage 
(Wiltshire et al. 2013) - is excluded. The data in the range 
6.25 < r ^ 12.5/i“^Mpc common to both the hrst primed 
and unprimed shells is important in establishing the boost 
offset which is more pronounced at small r. To account for 
systematic uncertainties, we have therefore applied a con¬ 
tinuous variation of the first shell boundary in the range 2 
- 6.35 h“^Mpc, while keeping the widths of the shells fixed. 
This leads to a value of p = —0.88 ± (0.25)stat ± (0.13)sys 
where the hrst uncertainty is the statistical and the second 
systematic. For the case of the primed shells we also note 
the corresponding velocity calculated from the best ht value 
of A is n = 772^440^ km sec“^, which is indeed close to the 
actual boost magnitude of 635 ± 38 km sec“^, albeit with a 
very large uncertainty. 

We repeated the analysis using 8 shells rather than 11 
to smooth out variations that could interfere with the sys¬ 
tematic boost ohset. The second conhguration uses shells 
of width 18.75 /i“^Mpc, starting from an inner cutoh of 
2 h“^Mpc and 9.375 h“^Mpc for unprimed and primed shells 
respectively. We hnd p = —0.89 ± 0.34 for the primed shells 
and p = —0.96 ± 0.26 for the unprimed shells. With a 
continuous variation of the inner shell boundary from 2 - 
9.375 h“^Mpc, we arrive at a value ofp = — 0.87±(0.33)stat± 
(0.09)sys. Fig. 2(b) shows the resultant best ht curves. 

In Fig. 2 we note a discrepancy between the best ht 
power law and the negative values on shells in the range 
of 40 - 60h“^Mpc (or (rf) = 1600-3600(/i“^Mpc)^). This 
is understood to be the result of structures in this particu¬ 
lar range, which give rise to both a residual monopole and 
dipole variation of the Hubble expansion in the LG frame. As 
Wiltshire et al. (2013) show, there is evidence that the boost 
to the CMB frame somewhat compensates for structures in 
this range. One hnds that in the range 40 < r < 60 h“^Mpc 
(and only in this range) the monopole variation is less in 
the CMB frame, {SHa)faj^^ < {SHa)^^^, while the dipole 
magnitude is also less in the CMB frame, becoming consis¬ 
tent with zero in the middle of the range. If the boost to the 
CMB frame exactly compensated for structures in the range 
40 < r < 60/i“^Mpc then the dipole magnitude should re¬ 
main close to zero in shells at larger distances. However, the 
magnitude of the CMB frame dipole increases to a resid¬ 
ual offset at large distances, while the LG frame dipole is 
consistent with zero in most outer shells. 

Thus it appears that the dipole almost - but not en- 
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Figure 2. Best fit power law to the radial variation in the spherically averaged Hubble law in 2 configurations of; (a) 11 shells; (b) 8 
shells. The dashed blue (lower) curve indicates the best fit to primed shells only (empty circles), the dotted red (upper) curve indicates 
the best fit to unprimed shells only (filled circles). The solid black curve is the combined best fit using both primed and unprimed shells. 
The first data point - corresponding to unprimed shell 1 - is omitted in each plot, as it is off the scale. 


tirely - has the character of a Lorentz boost dipole. Given 
that there are structures that give rise to a residual nonlin¬ 
ear Hubble expansion in the range 40 < r < 60h“^Mpc we 
cannot expect a perfect power law fit (15) in Fig. 2. How¬ 
ever, the deviation from a power law is consistent with the 
observation that in the range over which 

the boost almost compensates for nonlinear structures. 

Now that we have verified the power law nature of the 
difference we must check whether this re¬ 

sult is unique for the boost to the LG frame. To investi¬ 
gate this we determine the Hubble constant in radial shells 
for frames boosted arbitrarily with respect to the CMB, 
denoted by frame “X”, and then fit (15) to the resulting 
AH = ^CMB ~ curve. We vary the direction of the 
boost to frame X while holding the magnitude constant, 
thus producing a sky map. We first choose a magnitude of 
635 km sec“^ corresponding to the boost from the GMB to 
LG frame of reference. 

To display these sky maps in a meaningful fashion we 
cannot simply plot the value of p. Suppose that the GMB 
is boosted from a frame which has AHs = A{(r1 ')with 
p = — 1 and A > 0, representing the best fit boost offset. If 
one now boosts in the opposite direction by 635 km sec“^ 
then one finds a best fit power law with AHs = A((rf )^)^ 
with p « — 1 but A < 0 since the GMB frame necessarily 
has the smaller value of Hs on average. In each case we must 
first of all determine whether (15) gives a better overall fit 
with A>0orA<0 - given that some data points will 
always be opposite to the overall trend. In Fig. 3 we plot® 


r |p+l|, A 5=0 

\ 2-|p+l|, A<0 


(16) 


which takes the value /p = 0 at the best fit with A > 0 and 
fp — 2 at the best fit with A < 0. The latter point turns out 


® Both primed and unprimed shells are used (in the 11 shell case), 
to produce a smoothed sky map without determining systematic 
uncertainties. 


to be in the opposite direction, but not exactly opposite the 
best fit direction, reflecting the uncertainties in the method. 

Fig. 3 shows that for a boost magnitude of 635 km sec“^ 
the Local Group is indeed contained in a set of frames that 
display strong evidence of being a minimum Hubble expan¬ 
sion variation frame as p « —1. More precisely, the monopole 
Hubble expansion in the CMB frame compared to frames 
boosted at 635 km sec“^ relative to the CMB has the math¬ 
ematical character of a systematic boost offset for directions 
close to that of the LG. In Fig. 3 we can see a distinctive dif¬ 
ference between the directions for which the boosted frame 
has the lesser variation (A > 0, values plotted closest to 0), 
and the directions for which the CMB frame has the lesser 
variation (A < 0, values plotted closest to 2). 

We have verified that the additional monopole variation 
of the Hubble expansion seen in the CMB frame does have 
the character of what is expected by a boost from the LG 
frame, if the LG frame is close to the frame in which the 
monopole variation is minimized. However, as yet we have 
not varied the magnitude of the boost. Given the uncertain¬ 
ties we have already noted, it is clear that any systematic 
boost offset will become very hard to confirm statistically 
if the boost magnitude is small, since the other uncertain¬ 
ties will then become dominant. The systematic boost offset 
method can only give a rough indication of the boost di¬ 
rection for the large boosts rather than defining a precise 
“minimum Hubble expansion variation frame”. There are in 
fact many degenerate frames. 


4 MINIMIZING THE AVERAGE SPHERICAL 
HUBBLE EXPANSION VARIATION 

We will now investigate arbitrary variations of the boost 
magnitude and direction to determine to what extent we can 
define a frame of reference in which the monopole variation 
in the Hubble expansion is minimized. 
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0.0 0.3 0.6 0.9 1.2 1.5 1.8 

fp— Ip + 11 when ^ > 0 and 2 — |p + 11 when A < 0 

Figure 3. The best fit parameters to a systematic boost offset (15) for frames boosted from the CMB frame at 635 km sec“^. The black 
cross denotes the boost to the LG frame and the white cross denotes the boost to the frame with minimum variation in the spherically 
averaged Hubble law for a boost of this magnitude, which will be discussed in §4.1. In all figures, the galactic longitudes ^ = 0°, 180°, 360° 
are on the right edge, centre and left edge respectively. 


4.1 Variations in the “nonlinear regime” 


Initially we will consider monopole variations with respect 
to a uniform 5H = 0 expectation below the scale of sta¬ 
tistical homogeneity (< 100/i“^Mpc). This is quantified by 
summing the mean square differences of (9), to give a statis¬ 
tic 


xl{nf,ni) = 




Hi-fi+H?-% 


(17) 


where Uf and n; define the upper and lower shells included 
in the range of the calculation respectively. We will take the 
primed shells with, rii = 1' and n/ = 8', covering the range 
6.25 < r ^ 106.25/i“^Mpc. This includes all data poten¬ 
tially in the regime of nonlinear Hubble expansion, while 
excluding the innermost unprimed shell which may have in¬ 
complete sky coverage. 

In identifying a frame of “minimum nonlinear Hubble 
variation” we must also estimate an uncertainty about this 
minimum. Since we do not expect an exact fit to a linear 
Hubble expansion in any frame it is inappropriate to use the 
ordinary confidence intervals of the Xa distribution, as this 
would not give a fair measure of the improvement in unifor¬ 
mity between two reference frames, instead rendering all ref¬ 
erence frames at least® 1 a from the “expectation”. Instead 
we follow Wiltshire et al. (2013) in using Bayesian statistics 
to give a measure of the relative probability of uniformity 
of Hubble expansion in different reference frames. In partic¬ 
ular, using the complementary incomplete gamma function 
for the Xa distribution we directly calculate the probabilities 
Pjj, of the spherically averaged Hubble law giving a value 
that coincides with H^, using the shells 1' to 8', for the pair 
of reference frames U and V. A Bayes factor B = Py/Py 
then computed. We use the bounds InB = 1, 3, 5 to denote 


® At the global minimum we find a Xa P®r degree of freedom of 
2 . 0 . 


confidence regions, where the range 1 < InB ^ 3 represents 
positive Bayesian evidence, 3 < In B ^ 5 represents strong 
Bayesian evidence and In B > 5 represents very strong ev¬ 
idence (Kass & Raftery 1995). Thus if we take frame U to 
be the frame with the minimum Xa within a specified set 
of reference frames, any boosted frame V with Bayes factor 
InB ^ 1 relative to U can be considered statistically equiv¬ 
alent to U in the sense that the difference in probability 
is “not worth more than a bare mention” (Kass & Raftery 
1995). 

The frame of reference with the minimum monopole 
variation is found using a downhill optimization with (17). 
This reveals a global minimum variation (MV) frame for 
a boost in the direction {£,b) = (59.3°, 16.6°) with magni¬ 
tude 740.61487.0 km sec ^ with respect to the LG frame. 
(The uncertainty is taken from the InB ^ 1 bound.) Even 
though the boost is large, the uncertainties are also large. 
The Bayes factor In {Pmv/Plg) ~ constitutes positive 
but not strong evidence for the MV frame relative to the LG 
frame. However, the corresponding boost from the CMB ref¬ 
erence frame to the global minimum frame has a magnitude 
1203+111 km sec which is even larger and now very sig¬ 
nificantly nonzero, with In (^mv/^cmb) ~ 33.3, which is 
extremely strong evidence against the GMB frame. 

To visualize the distribution of fit in the 3-dimensional 
parameter space {v, i, 6} we show two angular slices at fixed 
values of v in Fig. 4, and a slice along the locus of {£, b) 
values for which Xa is minimized for fixed v in Fig. 5. The 
angular distribution in Fig. 4 is similar for both velocities. 
However, the Bayesian confidence regions grow large as v is 
decreased, eventually taking up the whole sky for very small 
velocities. Only for large boosts is there a well-defined boost 
direction that reduces monopole variation. 

In Fig. 5 for each boost magnitude we locate the direc¬ 
tion of minimum Xa and plot the corresponding value. The 
angular coordinates (£, b) of the minimum are different in 
each case. The distribution of Xa relative to the InB 1 
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Figure 4. Contour maps of angular variation of for two choices of boost magnitude with respect to the Local Group: (a) 740 km sec~t 
(b) 635km sec“t. The solid contours show the levels of Bayesian evidence InB = P{t.,b)) for each direction with respect 

to the frame of minimum Xa in each case. The white cross on (b) shows the direction of the boost to the CMB frame (also of magnitude 
^ 635km sec“t), Jn all figures, the galactic longitudes 1 = 0°, 180°,360° are on the right edge, centre and left edge respectively. 



Boost velocity (km 


Figure 5. Variation of the minimum Xa with fixed boost velocity, v. The locus of (£, b) values for which this minimum is found lies 
within Ri 20° of the galactic plane for cases in which InB = ln(P|^y/P^) 1 as indicated by the dashed line (where MV is the global 

minimum and V any frame along the locus). 


confidence interval (dashed vertical lines on Fig. 5(a)) re¬ 
veals the main problem in constraining a “minimum varia¬ 
tion frame” by this technique. The near flat distribution of 
Xa values within the In B $( 1 bound of the global minimum 
are found to lie on a locus of boost directions {i, b) which 
all lie close to the galactic plane. This is of course the Zone 
of Avoidance region, where the COMPOSITE sample lacks 
data, due to the Milky Way obscuring distant galaxies. Evi¬ 
dently, we are free to perform large boosts in the plane of the 
galaxy, as the data is not constrained there. This hypothesis 
could be checked by simulating data with the same charac¬ 
teristics as the COMPOSITE sample, using exact solutions 
of Einstein’s equations (Bolejko et al. 2015). Such an inves¬ 


tigation is beyond the scope of this paper. For now, with 
the available data we can only conclude that the LG frame 
is not ruled out as the standard of rest by this criterion, 
since the Bayesian distinction between the two frames with 
I lnB| = 2.77 is not strong. 

We have found a set of degenerate frames of reference 
which might be taken as the minimum average monopole 
variation frame. Although large uncertainties still exist, we 
are able to see that the boost to the CMB frame, shown by 
a white cross in Fig. 4(a), is far from our degenerate set of 
possible boosts to the minimum variation frame. In Table 1 
we make a comparison similar to Table 1 of Wiltshire et al. 
(2013) between the LG frame and the global minimum xi 
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Table 1. Hubble expansion variation in radial shells in minimum Hubble expansion variation (MV) and LG frames. Spherical averages 
(4) are computed for two different choices of shells, rs < r ^ ^s+li the second choice being labeled by primes. In each case we tabulate 
the inner shell radius, the weighted mean distance, fs; the shell Hubble constants, {Hs)^q and (Hs)j^jy in the LG and MV frames, 
and their uncertainties determined by linear regression within each shell, together with its “goodness of fit” probability Qs and reduced 
(for = Ng — 1); InB where B is the Bayes factor for the relative probability that the MV frame has more uniform SHg = 0 than 


frame when is 

summed in 

all shells with r > 

S • 

and cTs 

are given 

in units hkm sec 

^ Mpc 

■1 


Shell s 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

Ng 

92 

505 

514 

731 

819 

562 

414 

304 

222 

280 

91 

Tg ( h~^Mpc) 

2.00 

12.50 

25.00 

37.50 

50.00 

62.50 

75.00 

87.50 

100.00 

112.50 

156.25 

fg ( h^^Mpc) 

5.43 

16.33 

30.18 

44.48 

55.12 

69.24 

81.06 

93.75 

105.04 

126.27 

182.59 

{Hg)LG 

117.9 

103.1 

106.5 

105.5 

104.8 

102.1 

102.8 

103.2 

103.7 

102.4 

101.0 

Ds)lG 

4.6 

1.4 

1.0 

0.7 

0.7 

0.7 

0.9 

0.9 

1.0 

0.8 

1.7 

(Qs)lG 

0.000 

0.000 

0.000 

0.000 

0.998 

0.940 

1.000 

1.000 

1.000 

0.993 

0.999 

ixi/l'g) LG 

23.656 

7.767 

2.185 

1.419 

0.864 

0.909 

0.594 

0.542 

0.622 

0.803 

0.590 

{Hs)mV 

118.5 

102.9 

106.7 

104.5 

104.8 

102.9 

102.6 

103.9 

104.9 

102.7 

102.0 

(^s)mV 

4.6 

1.4 

1.1 

0.7 

0.7 

0.7 

0.9 

0.9 

1.0 

0.8 

1.7 

{.Qs)mv 

0.000 

0.000 

0.000 

0.000 

0.330 

0.887 

1.000 

1.000 

1.000 

0.964 

0.999 

ixllVg)MV 

29.130 

12.320 

3.037 

2.005 

1.021 

0.928 

0.682 

0.600 

0.667 

0.854 

0.603 

In B (r ^ Tg) 

3.53 

2.85 

2.79 

1.99 

0.85 

0.33 

0.36 

0.14 

0.07 

0.45 


Shell s 

1' 

2' 

3' 

4' 

5' 

6' 

7' 

8' 

9' 

10' 

11 

Ng 

321 

513 

553 

893 

681 

485 

343 

273 

164 

206 

91 

Ds (/l“^Mpc) 

6.25 

18.75 

31.25 

43.75 

56.25 

68.75 

81.25 

93.75 

106.25 

118.75 

156.25 

(h“^Mpc) 

12.26 

23.46 

37.61 

49.11 

61.74 

73.92 

87.15 

99.12 

111.95 

131.49 

182.59 

(Hs)lG 

103.5 

103.5 

103.9 

106.6 

103.9 

102.0 

103.2 

103.6 

101.6 

102.7 

101.0 

(CTs)iG 

1.8 

1.1 

0.9 

0.7 

0.8 

0.8 

0.9 

0.9 

1.0 

0.9 

1.7 

{.Qs)lG 

0.000 

0.000 

0.000 

0.031 

0.960 

1.000 

1.000 

1.000 

0.996 

0.999 

0.999 


11.427 

3.246 

1.792 

1.090 

0.907 

0.701 

0.592 

0.608 

0.728 

0.711 

0.590 

{Hg)MV 

102.7 

104.3 

103.3 

106.1 

104.2 

102.9 

102.9 

104.7 

102.7 

102.9 

102.0 

(cts)mV 

1.8 

1.1 

0.9 

0.7 

0.8 

0.8 

0.9 

0.9 

1.0 

0.9 

1.7 

{Qs)mV 

0.000 

0.000 

0.000 

0.000 

0.481 

1.000 

1.000 

1.000 

0.967 

0.997 

0.999 

{xl/^g)MV 

18.547 

4.940 

2.429 

1.428 

1.002 

0.734 

0.704 

0.613 

0.807 

0.752 

0.603 

InB (r ^ Tg) 

2.38 

2.32 

2.34 

1.93 

0.55 

0.33 

0.41 

0.14 

0.24 

0.50 



(MV) frame. In particular, we directly calculate the prob¬ 
abilities Pj^Y spherically averaged Hubble law 

giving a value that coincides with , when all shells larger 
than a given shell are included. If we exclude the inner¬ 
most unprimed shell 1, then the Bayesian evidence for the 
MV frame having a more uniform average Hubble expansion 
than the LG frame, lies at most in the range 1 < InH $( 3, 
which represents positive but not strong Bayesian evidence 
(Kass & Raftery 1995). By contrast, the Bayesian evidence 
that the LG frame has a more uniform spherically average 
Hubble expansion than the CMB frame is very strong with 
InH > 5 (Wiltshire et al. 2013). 


4.2 Variations in the “linear regime” 

The goodness of fit, Qg, of a linear Hubble law in the in¬ 
nermost shells of Table 1 is poor, as we would expect since 
these shells are in the nonlinear regime. Beyond approxi¬ 
mately 75/i“^Mpc we expect to pass into the linear regime 


(Scrimgeour et al. 2012), and this is seen with the decreasing 
values of fifjvg, where Xg is given by (3) for each shell. 

However, contrary to expectation, the MV frame deter¬ 
mined from (17) has an overall poorer goodness of fit Qg 
than the LG frame to a linear Hubble law in shells with 
s 4', including in particular in shells 9' and 10 which 
should be in the linear regime. We note that the asymptotic 
value, Hq, is 1% higher in the MV frame as compared to the 
LG frame, and this may contribute to SHg being smaller, 
even though the goodness of fit in individual shells is poorer 
in some cases^. 

Any true candidate for a minimum Hubble expansion 
variation frame should also clearly demonstrate the emer¬ 
gence of a linear Hubble law consistent with the existence of 
a statistical homogeneity scale. The Xa statistic (17) involves 

^ The value of is also 1% larger in the LG frame as compared 
to the CMB frame, but by contrast the goodness of fit in the 
linear regime shells is better in the LG frame than in the CMB 
frame (Wiltshire et al. 2013). 
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Figure 6. Variation of the minimum Xb/^ with fixed boost velocity, v. 



(a) 



(b) 


Figure 7. (a) Variation of Xa along a locus of (£, b) values for which Xb minimized at each fixed boost velocity v. (b) Variation of 
Xb/^b along a locus of (£, b) values for which Xa minimized at each fixed boost velocity v. 


minimizing the variation SHs = {Hs — relative to 

the asymptotic Hubble constant. But a boost can also alter 
Up in a way which makes for a worse goodness of fit to a 
linear Hubble law. Therefore it is not a completely suitable 
candidate statistic. 

In order to quantify the emergence of a linear Hubble 
law we will therefore alternatively minimize the quantity 

xl = ^xl ( 18 ) 

where Xs is given by (3) in the sth shell. This sum is per¬ 
formed over the unprimed configuration of shells as shell 
7 has an inner cutoff near the boundary of the nonlinear 
and linear regimes. Thus (18) gives a measure of the good¬ 
ness of fit to a linear Hubble law averaged over the outer 5 
shells, without normalizing the asymptotic Hubble constant. 
Defining the total degrees of freedom Uh = ~ 1> 

we find Xb/^b = 0.631 for the LG frame, XhKh = 0.692 
for the MV (minimum Xa) frame, and Xb/i'b = 0.653 for 


the CMB frame. Thus even the CMB frame shows a clearer 
emerging linear Hubble law than the minimum Xa frame. 

To determine whether there is any frame with a more 
definitive emerging linear Hubble law than the LG frame, 
we determine the distribution of (18) upon making arbitrary 
boosts with respect to the LG frame. 

In Fig. 6 we locate the direction of minimum Xb af each 
boost magnitude and plot the corresponding value, analo¬ 
gously to Fig. 5. We find a best fit boost of 222.3 km sec“^ 
in the direction {£, b) = (241.84°, 70.53°) with respect to the 
LG frame, with a value of Xb/*^b = 0.621. This is 45° from 
the direction of the residual CMB temperature dipole in the 
LG frame, and so does not appear related. 

In Fig. 7 for each boost magnitude we calculate the 
value of Xa in an (£, b) direction determined by minimizing 
with respect to Xm nnd vice versa. Thus, we compute the 
locus of {l,b) values in the {v,l,b} parameter space that 
minimize Xb for each fixed v, and then compute Xa at these 
parameter values, and vice versa. It is apparent that making 
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boosts of the order > 100 km sec ^ along the locus of {£, b) 
values which minimize Xh results in an increase in Xa to 
values which are disfavoured by Bayesian evidence. This is 
because making improvements in Xb requires boosts in {£, b) 
directions away from our minimum Xa- 

However, we note that: (i) all values of Xb shown in 
Fig. 6 and 7 are consistent with very probable fits, given the 
low values of Xb/^u; (ii) there are less data points per se in 
those outer shells used in Xb s-s compared to the inner shells 
used in Xai (iii) local boosts have a relatively small 
effect on the values of czi in the outer shell. These facts 
mean we should also exercise caution about drawing strong 
conclusion solely from the minimization of Xb- The primary 
utility of the Xb statistic is not as a sole discriminator, but 
as a check against the potential bias in the Xa statistic due 
to its anchoring to the value in the outermost shell. 

Finding a joint minimum for Xa a-nd Xb would be feasi¬ 
ble if the {£, b) values of each global minimum were close on 
the sky. However, this is not the case - the global minimum 
for Xa is at an angle of 93° from the global minimum of Xm 
i.e., they are roughly orthogonal. As there is no unbiased 
way to weight these two statistics we cannot set out to de¬ 
termine a weighted minimum as any result would be highly 
sensitive to our choice of weighting. 

4.3 Systematic boost offsets from the Local Group 

Neither statistic Xa nor Xb appears entirely satisfactory for 
establishing a global minimum expansion variation frame. 
The Xa statistic is the better measure of Hubble expansion 
variation in the nonlinear regime but is also affected by po¬ 
tential bias in the anchoring of H^. The most we can say 
is that there is a freedom to perform large boosts in the 
plane of the galaxy, given the lack of data in the Zone of 
Avoidance. 

If Xa is taken as the better statistic, then a criterion for 
breaking the boost degeneracy may be possible by returning 
to the systematic boost offset analysis of §3. Any true best fit 
frame should show a clear signal of a boost offset (12) with 
respect to the Local Group. The “best” boost offset can be 
characterized in 3 ways, each with its own challenges: 

(1) Determine the boost for which p = —1. This is hindered by 
the fact that there are many boosts that satisfy this crite¬ 
rion, at almost every magnitude from the LG. 

(2) From the value of A in (15) determine a derived boost ve¬ 
locity, Uder. Any boost offset should have Uder consistent with 
the true boost magnitude Vtrue within uncertainties. How¬ 
ever, this is difficult due to the large uncertainties associated 
with the value of A. 

(3) Determine a measure of variation in the fit of the boost 
offset, given by (AS) in the Appendix. This is also prob¬ 
lematic since all fits are extremely good due to the large 
uncertainties in the Ha. 

We will therefore use method (1) to determine the di¬ 
rection of the boost on the sky, and then consider (2) and 
(3) to constrain the magnitude®. 

® The choice of shell boundaries introduces systematic uncertain¬ 
ties in this analysis. In the production of the sky maps in Fig. 8 
we calculate both primed and unprimed shells and fit the power 


First we check for a systematic boost offset for the 
global Xa minimum frame determined in §4.1. A sky map 
of fp values as given by (16) for boosts of magnitude 
740km sec“^ is given in Fig. 8(a), with the InB ^ 1 and 
1 < InH ^ 3 confidence regions for Xa displayed. We note 
that there are in fact boosts with values of p « —1 consis¬ 
tent with the Xa minimum. However, these directions are far 
more constrained and do not align with the exact minimum. 
In addition, the {£, b) direction at this magnitude with the 
best fit to (15) has an inconsistent value of the derived ve¬ 
locity. Thus we do not see a clear systematic boost offset 
between the Local Group and the frame corresponding to 
the global minimum xL further ruling this out as a poten¬ 
tial candidate for the standard of rest we are looking for. 

The next step is a global search for the best systematic 
boost offset from the LG. In order to further understand 
the angular distribution of fp values (16) for boosts from 
the LG frame, we arbitrarily choose a boost magnitude of 
200km sec“^ and plot fp with respect to {£,b) in Fig. 8(b). 
We have found that for all interpolating velocities between 
the 200 km sec“^ and 740 km sec“^ cases displayed in Fig¬ 
ure 8 there is a region of {£, b) values for which fp « 0. Thus 
in order to use this method to find a realistic systematic 
boost offset we must use an additional criterion. 

In Fig. 9 we plot the values of Uder/utrue and S/{n — 2) 
(for n data points) from (A8) with respect to the boost mag¬ 
nitude, where for each magnitude the {£, b) direction is that 
for which p is closest to —1 and A > 0 (i.e. fp « 0). Thus, 
we can use these additional quantities to constrain a sys¬ 
tematic boost offset along the locus of {£, b) directions in the 
3-dimensional {u, £, b} parameter space. The expected value 
of S has a x^ distribution for (n — 2) degrees of freedom, and 
thus S/{n — 2) has an expectation value of unity (York et al. 
2004). Glearly, the values of S/{n — 2) in Fig. 9(b) are consis¬ 
tent with a very good fit to (15) for all boosts. Our inability 
to tightly constrain the boost magnitude is no doubt due to 
the lack of data in the Zone of Avoidance and large uncer¬ 
tainties in the values of Although (A8) it is not 

useful for constraining the boost magnitude, we nonetheless 
see that the ratio of derived and true velocities in Fig. 9 does 
show a meaningful difference on this interval. 

Using both primed and unprimed shells in our calcula¬ 
tion we find that the Uder/utrue = 1 at Utrue = 171km sec“^ 
in a direction {£,b) = (59°,— 4°) ± (5°, 4°). For this boost, 
we find the result p = —1.18 ± (0.95)stat ± (0.32)sys, where 
the systematic uncertainty is determined as in §3, consis¬ 
tent with p — —1. The difference between the spherically 
averaged Hubble expansion in this frame, which we denote 
by X and the LG is shown in Fig. 10. One can see that a 
systematic boost offset is apparent here. However, due to 
the large uncertainties that arise when taking differences of 
the Ha the result is also statistically consistent with zero. 
(Thus the question of whether the first unprimed shell 1 
should be included in the analysis, due to its incomplete sky 
cover, is immaterial. Indeed, if we take only primed shells 
then we obtain p = —1.17 ± (1.2)atat.) This frame X is only 
weakly disfavoured compared to the global Xa minimum. 


law (15) to all points. For Fig. 9 we consider fits to primed and 
unprimed only, and to both. 
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o.'o 0^3 0^6 0^9 1.2 ' i!5 L8 2 !i 

fp~ Ip + 11 when A>Q and 2 — |p + 11 when A < 0 


Figure 8. The best fit parameters to a systematic boost offset for boosts from the Local Group of magnitude: (a) 740 km sec“^; (b) 
200 km sec“^. The thick blue contours denote the corresponding Xa distribution. The solid contours show the levels of Bayesian evidence 
\nB = b)) “ each direction (^, 6) with respect to the frame of minimum Xa each boost magnitude. In (a) 

In B = 1, 3, 5 contours are visible while in (b) only the In B = 1, 3 contours are visible, as no directions have In B > 5 in that case. In all 
figures, the galactic longitudes ^ = 0°, 180°, 360° are on the right edge, centre and left edge respectively. 




Boost velocity (km s 0 


(b) 


Figure 9. (a) The ratio of the derived velocity from the best fit power law and the true boost velocity, at each v the direction is 
determined by the best fit for which p is closest to —1. (b) The statistical variation in the same direction as determined by S/(n — 2) 
from (A8), where n is the number of data points being fit. 


ln(Pj^y/P,(-) = 1.7, and is within Ict of the global minimum 
of Xbi similar to the LG frame on both counts. 

Our choice of frame X above is based on taking 
"Uder/ftrue = 1, a couditiou which may only be approximately 
matched in reality, given our huge uncertainties in the co¬ 
efficient A. We again have a degeneracy in the choice of 
minimum Hubble expansion variation frame that satisfies 
the two conditions p = —1 and Uder = Utrue. 

The results of this section confirm the finding of §4.1 
that our determination of a suitable cosmic rest frame is 
limited in the COMPOSITE sample by a degeneracy un¬ 


der boosts close to the plane of the galaxy. The consistency 
between the methods of this section and §4.1 may be less 
significant, as they are not completely independent. In the 
LG frame the primary source of the monopole variation is 
the increased value of Ha in the innermost shells, while the 
more distant shells closer to the linear regime show closer to 
asymptotic values. Thus boosting to a frame with a reduced 
Ha in the innermost shells will give the most significant im¬ 
provement to xi, relative to which small changes in the more 
distant shells are negligible. This is precisely the type of dif¬ 
ference we model with a power law of the form (15) with 
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Figure 10. The variation (9) in the spherically averaged Hubble law: (a) 5Hs in the LG frame; (b) SHs in frame X; (c) the systematic 
boost offset between the LG and frame X. (Note: the first unprimed shell 1 is shown.) 


p « — 1. Consequently, if our hypothesis concerning (15) is 
correct then it is not surprising that we see the consistency 
in the angular directions that minimize on one hand, and 
which give values of p « —1 with A > 0 on the other. 


4.4 Angular Hubble expansion variation 

Wiltshire et al. (2013) also explore the extent to which an¬ 
gular averages of the Hubble expansion offer an independent 
characterization of a minimum Hubble variation frame of 
reference. When one takes angular Gaussian window aver¬ 
ages of the Hubble expansion a dipole becomes apparent. 
Wiltshire et al. (2013) show that this dipole is strongly cor¬ 
related with the residual CMB temperature dipole when 
both are referred to the LG (or LS) rest frame. If we 
are to define a new cosmic standard of rest, within which 
we still observe a residual CMB temperature dipole, then 
Wiltshire et al. (2013) argue that such a dipole must have 
a nonkinematic origin. The correlation of the residual CMB 
temperature dipole and Hubble expansion dipole supports 
the proposal that structures in the nonlinear regime of 


expansion are simultaneously responsible for these effects. 
Thus we are interested in finding a frame of reference in 
which this correlation is maximized. 

We have investigated this question, and find that cor¬ 
relation of the residual CMB temperature and Hubble vari¬ 
ation dipoles under arbitrary boosts does not offer a viable 
characterization of the minimum variation frame (McKay 
2015). By making boosts in the appropriate direction from 
the LG we are able to artificially increase both the magni¬ 
tude of the Hubble expansion dipole and the CMB residual 
temperature dipole simultaneously, at the expense of also in¬ 
creasing the monopole expansion variation. By calculating 
the Hubble expansion dipole and higher multipole coeffi¬ 
cients using HEALPIX® in the frame of maximum correla¬ 
tion for a given boost magnitude from the LG we can ob¬ 
serve this artificial increase in the strength of the dipole rel¬ 
ative to the higher multipoles. Boosting in directions which 
make the dipoles more pronounced will naturally increase 
the correlation, but this is irrelevant if the monopole varia- 

® http://healpix.jpl. nasa.gov/(Gorski et al. 2005) 
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tion is also increased. Therefore, as this artificially induced 
increase in the correlation cannot be distinguished from a 
physically meaningful increase, this method does not offer a 
viable characterization of the minimum variation frame and 
was abandoned as a line of investigation (McKay 2015). 


5 HUBBLE EXPANSION VARIATION IN THE 

COSMICFLOWS-2 CATALOGUE 

Thus far our investigation has been based entirely on the 
COMPOSITE catalogue of distances and redshifts. In this 
section we aim to repeat the monopole Hubble expansion 
variation analysis on the recently released Cosmicflows-2 
(CF2) catalogue. Systematic differences become apparent in 
this analysis which we will investigate in §6. 

CF2 is a compilation of distances and redshifts from 
both new and existing sources of observational data. The 
entire catalogue consists of over 8000 galaxies both locally 
and extending beyond the scale of statistical homogeneity. 
This compilation has a large subset of galaxies and galaxy 
clusters in common with the COMPOSITE sample, includ¬ 
ing the large SFlH—h sample of Springob et al. (2007). In 
total, the distances are determined by six different methods 
and compiled together by Tully et al. (2013). The CF2 data 
are presented in two sets, one with all iirdividual galaxies 
included, and one condensed into galaxy groups, including 
groups consisting of one galaxy. We will use the entire data 
set of 8162 galaxy redshifts and distances, freely available 
from the extragalactic distance database^®. 

5.1 Spherically averaged Hubble expansion and 
treatment of biases 

We repeat our earlier analysis to calculate SHs (9) in the 
same two configurations of 11 spherical shells. The CF2 data 
are presented with a modified “recession velocity”, cz, and 
a raw observed redshift. Given the prevalence of the use 
of such modifications, particularly in bulk flow studies, it 
is worthwhile to briefly investigate the effect such modifi¬ 
cations can have on monopole Hubble expansion variation. 
Tully et al. (2013) define an adjustment assuming a FLRW 
cosmology with = 0.27 and = 0.73. This adjust¬ 
ment is given by a Taylor expansion to 0{z^) of a homoge¬ 
neous isotropic expansion law, 

czmod = cz [l -I- i(l - qo)z - i(l - (jp - 3(7^ -I- l)z^] , (19) 

where — 20^^) and z is the redshift in the 

CMB frame. 

Since we wish to deal with cosmological model- 
independent quantities, this is not the type of adjustment 
that should be made. In particular, a homogeneous isotropic 
expansion law cannot be assumed below the scale of statisti¬ 
cal homogeneity if the conclusions of Wiltshire et al. (2013) 
are correct. Nor should such an expansion law be assumed in 
the CMB rest frame. However, for completeness we consider 
the adjustment (19) in order to rule it out as the cause of 
much larger systematic differences we will discuss shortly. 

Cosmicflows-2 distances retrieved 14 October 2014 from 
http://edd.ifa.hawaii.edu/. 


Fig. 11 shows SH using the CF2 sample. It is imme¬ 
diately evident that these plots are very different to those 
found with the COMPOSITE sample, as given in Fig. 3 
of Wiltshire et al. (2013). Most notably, the uniformity of 
Hubble expansion in spherical shells is considerably worse 
in CF2 as compared to the COMPOSITE sample, in both 
rest frames. While there is a small shift introduced by the 
redshift modification, it is only makes a modest impact on 
the very large differences from a uniform expansion seen in 
the unadjusted redshifts. 

It should be noted that by (9) SHs depends on the 
asymptotically normalized Hubble constant, Hq. Thus the 
vertical shift seen in Fig. 11 when the FLRW “correction” 
is applied is primarily due to a change in the asymptotic 
value, Hq, as data at smaller values of cz are barely affected 
by the correction in (19). Table 2 gives the numerical val¬ 
ues of the Hubble constant with adjusted and raw redshifts, 
along with the number of objects in each shell and the mean 
shell radii, making it transparent how the shift in Figure 11 
arises. 

As we will see in §6, the treatment of Malmquist biases 
is the most likely cause of the systematic difference we see in 
Figure 11 as compared to Fig. 3 of Wiltshire et al. (2013) for 
the COMPOSITE sample. To study the nature of this sys¬ 
tematic difference we make use of the SFI-f-1- sample, which 
is a subset of both the COMPOSITE and CF2 samples. 
SFI-f-f (Springob et al. 2007) consists of Tully-Fisher (1977) 
Relation (TFR) derived distances for 4861 field and clus¬ 
ter galaxies. Since Springob et al. (2007) supply the SFI-f-|- 
sample with and without corrections for Malmquist biases, 
it makes it an ideal candidate with which to study the ef¬ 
fects that the treatment of such biases has on the spherically 
averaged Hubble expansion. 

To understand the effect of the Malmquist bias cor¬ 
rections applied by Springob et al. (2007), we calculate the 
monopole variation of the Hubble expansion for the SFI-I--I- 
sample with and without corrections. Figure 12 shows the 
significant difference in SHs between these two treatments. 

Since uncorrected SFI-I—I- data points are included in 
the CF2 catalogue, we can determine whether there is any 
systematic difference between this subsample and the re¬ 
mainder of the CF2 catalogue. If we take the CF2 catalogue 
and remove the 3625 points in common with our SFI-f-f 
sample, and repeat the analysis, we arrive at Fig. 13. We find 
that SHs does not change to any statistically significant ex¬ 
tent, by removal of this potentially biased data. This is an 
indication that the systematic bias present in the SFI-I--I- 
raw distances - uncorrected for Malmquist bias - is likely 
to also be present in the rest of the CF2 data. Thus we are 
confident that the discrepancy seen between the monopole 
variation of Fig. 11 for CF2 and Fig. 3 of Wiltshire et al. 
(2013) for the COMPOSITE sample is a systematic issue, 
arising from the treatment of Malmquist bias in the CF2 
catalogue, as we discuss further in §6. 

5.2 The systematic boost offset revisited 

One may ask whether, despite the obvious problems with a 
systematic bias in the CF2 data, any of the analyses applied 
to the COMPOSITE sample can nonetheless yield meaning¬ 
ful results. 

We applied the analysis of §4.1 in the CF2 catalogue. 
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Figure 11. The monopole Hubble expansion variation for the CF2 sample without the FLRW “correction” (19) (black filled circles) and 
with the FLRW “correction” (19) (blue crosses) in the: (a) Local Group frame of reference; (b) CMB frame of reference. 


Table 2. Hubble expansion variation in radial shells in CMB and LG frames for the CF2 data. Spherical averages (4) are computed for 
two different choices of shells, Vs < r ^ second choice being labeled by primes. In each case we tabulate the number of data 

points per shell, the weighted mean distance, fg; the shell Hubble constants, (Hs)i^q and (^s)cmb their associated uncertainties in 
the LG and CMB frames for both the raw redshifts and those adjusted with (19). 
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Figure 12. The monopole Hubble expansion variation for the SFI++ sample without corrections for Malmquist bias {black filled circles) 
and with corrections {blue crosses) in: (a) the Local Group frame of reference; (b) the CMB frame of reference. 
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Figure 13. The monopole Hubble expansion variation in the LG frame for: (a) all CF2 data {blue filled circles) and CF2 data without 
SFIH—h data {red crosses); (b) CF2 data without SFIH—h data {red crosses) and SFIH —\- only without Malmquist corrections {black filled 
circles). For the uncorrected SFIH—h data the first unprimed and primed values are not shown, these are 2.7 and 0.65 respectively. 


but found no reference frame in which Xa approaches unity, 
or even within the same order of magnitude. While some 
decrease in the variation of Hubble expansion was found 
for boosts in the galactic plane, the uncertainties were too 
large to give any statistically significant results. Such inves¬ 
tigations must be abandoned until the bias problems in the 
CF2 catalogue are dealt with. 

By contrast, we found that in spite of the bias problem, 
the signature of a systematic boost offset studied in §3 is 
nonetheless evident in CF2, as this involves the difference of 
the Hs values in the LG and CMB frames from Table 2, as 
plotted in panels (a) and (b) of Fig. 11. 

Fig. 14 shows the results of repeating the analysis used 
in §3, using the unadjusted CF2 distances. The best fit 
value for p in (15) is found to be p = —0.83 ± 0.17 for 
unprimed and p = —0.86 ± 0.26 for primed shells. Vary¬ 
ing the shell boundaries as in §3 we find a value of p = 
—0.84±(0.21)stat ±(0.06)sys. However, if we compare Fig. 14 
with Fig. 2, we see that there are more data points with 


which do not conform to the power 
law (15). However, the range of distances of the shells for 
which this is true coincides in Fig. 2 and Figure 14, being 
40< r < 60/i“^Mpc in the COMPOSITE sample and 
30 h~^ <r< 67/i“^Mpc in the CF2 sample. This is con¬ 
sistent with the hypothesis that aside from the systematic 
boost offset, there are structures responsible for nonlinear 
deviations in the monopole Hubble expansion in the range 
identified in the COMPOSITE sample, but untreated biases 
in the CF2 catalogue have broadened the range of distances 
attributed to the same structures'^. 

The systematic boost offset is still evident in the in¬ 
nermost shells of the CF2 sample. However, the fit is some¬ 
what worse than in the COMPOSITE sample due to more 


These data points are necessarily disregarded when we perform 
the logarithmic transformation to fit the power law (15), and so 
do not contribute to the stated uncertainties. 
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< rj > (h~^Mpc) 

Figure 14. Best fit power law to the radial variation in the spherically averaged Hubble law in 11 shells for the CF2 galaxies data. The 
dashed blue curve is the best fit to primed shells only (empty circles), the dotted red curve is the best fit to unprimed shells only (filled 
circles) and the solid black curve is the best fit to all data points. The first data point - corresponding to unprimed shell 1 - is omitted, 
as it is off the scale. 



U.O 0.3 0.6 0.9 1.2 1.5 1.8 

fp~ b + 11 when j 4 > 0 and 2 — |p + 11 when A < 0 

Figure 15. Best fit power law parameters to (15) across entire sky for boosts of 635 km sec“^ from the CMB frame. The cross indicates 
the direction of the boost to the LG, which is also of magnitude 635 km sec~^. In all figures, the galactic longitudes £ = 0°, 180°, 360° 
are on the right edge, centre and left edge respectively. 


data points lying in the increased range which deviates from 

(15) . Nonetheless, we can still check if the boost offset signa¬ 
ture is unique to the angular direction of the residual CMB 
dipole in the LG frame. Fig. 15 shows the value of fp from 

(16) which represent the best fit parameters for a system¬ 
atic boost offset for boosts of 635 km sec“^ across the en¬ 
tire sky. The results are consistent with those found for the 
COMPOSITE sample, providing further evidence that this 
is indeed not a random statistical outcome but is consistent 
with our hypothesis. 

We have also repeated the analysis of §4.3 for the best 
fit to a systematic boost offset from the LG frame, as shown 
in Figure 16. The angular directions found are consistent 
with our results for the COMPOSITE sample. For exam¬ 
ple, on the 200 km sec“^ sky map the best fit value is in 
the direction {l,b) — (55°, —5°), close to that found earlier. 
However, the value of p = —0.92 ±0.75 has a far greater sta¬ 


tistical uncertainty, which may well be due to the untreated 
biases. 

In conclusion, the greater number of data points in the 
CF2 catalogue may potentially yield statistically more ac¬ 
curate results than the COMPOSITE sample. However, at 
present this is prevented on account of untreated biases, to 
which we now turn. 


6 COSMICFLOWS-2 MALMQUIST 
TREATMENT 

The Cosmicflows-2 (CF2) and COMPOSITE catalogues 
deal with Malmquist distance biases in different ways. The 
results of the last section indicate that the treatment of such 
biases is crucial in establishing the actual nature of the vari¬ 
ation of cosmic expansion below the statistical homogeneity 
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fp— Ip + 11 when A>Q and 2 — |p + 11 when A < 0 


Figure 16. Best fit parameters to (15) using unadjusted CF2 data for boosts from the LG of: (a) 450 km sec (b) 200 km sec In 
all figures, the galactic longitudes £ = 0°, 180°, 360° are on the right edge, centre and left edge respectively. 


scale. Therefore we will perform further analyses to better 
understand these differences. 

Since the treatment of Malmquist bias is complex, we 
will first briefly remind the reader that in its current usage 
this term refers to at least three distinct biases that affect 
the average derived distance of a galaxy cluster: 

(i) Selection bias This is the simple homogeneous systematic 
error that results from using objects in a magnitude limited 
sample (Malmquist 1922). If the galaxies in a cluster have a 
true average flux (T’)p, since we can only observe the bright¬ 
est of these, we measure a biased average flux lF)b > 
with the effect of bias increasing with distance. Thus the 
luminosity distance dK = (L)^ /{Atv (F)) will be less when 
using a biased average flux, where (L)q is the average lumi¬ 
nosity as determined from standard candles calibrated using 
nearby objects. 

(ii) Homogeneous distribution bias This is a systematic av¬ 
erage error for objects around the same derived distance 
which can be understood in terms of statistical scatter. If 
we assume a standard Gaussian scatter, a, in the derived 
distances about an estimated mean, then - since the radial 
number density grows as N(r) oc r® - there are more objects 
with true distances larger than the estimated distance, than 
smaller. Thus at a given derived distance, more galaxies will 
have been scattered by the errors down from larger true dis¬ 
tances than up from smaller ones (Lynden-Bell et al. 1988; 
Hanski 1999). This is equivalent to giving more statistical 
weight to more distant values and thus the probability dis¬ 
tribution for the true distance is no longer Gaussian along 
the line of sight, centred on the measured distance, instead 
being skewed towards greater distances. For the special case 
of a constant density distribution, for example, one arrives 
at Eddington’s formula (Eddington 1914) 

E/(/rtrue l/Tder) — /^derived T 1.382(7 , (20) 

where p = 51ogr + 25 is the standard distance modulus 
when r is given in Mpc. 


(iii) Inhomogeneous distribution bias The inhomogeneous bias 
is analogous to the homogeneous one in that it involves a sys¬ 
tematic error in the statistical scatter of objects due to their 
distribution in 3-dimensional space. However, in this case it 
arises as number counts are higher in regions of greater den¬ 
sity, resulting in a systematic scatter of measurements out of 
higher into lower density regions (Strauss & Willick 1995). 
Thus the inhomogeneous effect is very sensitive to large vari¬ 
ations in large-scale structure along the line of sight. Failure 
to account for this type of bias can give spurious infall sig¬ 
natures on to high density regions. This bias is of course 
far more difficult to account for, requiring accurate density 
fields for structure along the line of sight for each observa¬ 
tion. 

The selection and homogeneous distribution biases typ¬ 
ically lead to underestimates of distances that increase as 
the distance grows, so that a plot of ptrue — Mder versus red- 
shift has a positive slope. However, the inhomogeneous dis¬ 
tribution bias can lead to the opposite effect. For example, 
Feast (1987) shows that in applying the TFR method when 
the spatial density of objects at a given 21cm line width is 
constant, then the required Malmquist correction is the clas¬ 
sical one given by Eddington (1914). However, when this is 
not the case it is possible to obtain overestimated distances 
(Feast 1987). 

Moreover, regardless of the details of any inhomoge¬ 
neous matter distribution, one just needs the spherically av¬ 
eraged N{r) to decrease sufficiently quickly for the direction 
of scatter in the standard homogeneous distribution bias to 
be reversed, giving overestimated distances. 

Some biases can be dealt with in the data reduction. In 
particular, by applying the inverse TFR method rather than 
the direct method one can in principle effectively eliminate 
the selection bias (Schechter 1980), leaving only a consider¬ 
ably smaller bias (Willick 1994, 1995). 

In the CF2 catalogue Tully et al. (2013) use an inverse 
TFR procedure to reduce the selection bias only, stating 
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that only a small subsequent correction for residual bias is 
required. In particular, they “make no adjustments for the 
distribution Malmquist effects” in their reported CF2 dis¬ 
tances (Tully et al. 2013). Their calibration carried out for 
this relation follows the procedures of Tully & Pierce (2000), 
Courtois & Tully (2012) and Sorce et al. (2013). 

On the other hand, the SFlH—h catalogue 
(Springob et al. 2007) which forms the major part of 
the COMPOSITE sample includes corrections to account 
for all homogeneous and inhomogeneous biases in their 
data. In their view their own treatment of the homogeneous 
and inhomogeneous distribution biases was “straightfor¬ 
ward” , given access to a reconstruction of the local density 
field (Erdogdu et al. 2006). However, they stated that their 
treatment of the selection bias was ad hoc because the 
selection criteria used are designed to mimic the observa¬ 
tional properties of the survey as closely as possible, and 
so are very inhomogeneous. They therefore provided both 
the raw and corrected distances, should other researchers 
adopt alternative methods for dealing with selection bias. 

The COMPOSITE sample incorporates the SFI-f-f 
distances with the Malmquist bias corrections of 
Springob et al. (2007), whereas CF2 uses a subset of 
uncorrected SFI-f-|- distances. While we do not have the 
data to independently repeat any of the bias corrections, 
we are able to test the consistency of assumptions made by 
Tully et al. (2013). 

6.1 The SFI-I—|- subsamples of CF2 

Tully et al. (2013) find that for 2071 common points be¬ 
tween their own survey and the SFI+-I- survey (excluding 5 
points judged to be “bad”) there was a “correction” of the 
form 

Afii = 0.492(±0.011) -h 0.000031(±0.000002 )c2ls (21) 

where A/ri = ficf 2 — Usfi, Ucf 2 is the CF2 distance mod¬ 
ulus with the zero point established by Courtois & Tully 
(2012), and fisfi is the Springob et al. (2007) unadjusted 
modulus with a nominal zero point consistent with Hq = 
100/ikm sec“^Mpc“^, and is the raw^^ redshift in the 
rest frame of the Local Sheet, which is close to the Local 
Group frame (Tully et al. 2008). We independently confirm 
the slope in (21) using the appropriate zero point^^, and plot 
this in Figure 17. Note that this comparison is made for a 
subset of the SFlH—h sample, henceforth SFI++A, consist¬ 
ing only of objects that are common between the SFI-I—|- 
survey and the independently obtained CF2 distances. In 
their final analysis Tully et al. (2013) use averages of the 
CF2 and SFI-I—I- distances with double weight given to the 
CF2 distances. 

The intercept in Fig. 17(a) simply reflects the different 
normalizations of Hq in the two samples. The positive slope 
of the linear relation in Fig. 17(a) is consistent with the CF2 
inverse TFR distances having a correction that accounts for 
the raw SFI-I—I- distances being increasingly underestimated 
due to the Malmquist selection bias. 

Let us now compare (17) to the Malmquist correction 

The FLRW adjustment (19) is not applied. 

The intercept in (21) is determined by a scaling of the data 


used by Springob et al. (2007) in the SFI-f-1- sample, and 
subsequently adopted in the COMPOSITE sample. We re¬ 
peat the analysis of Fig. 17 on the SFI-I—|-A subsample 
but now using the distances as corrected by Springob et al. 
(2007), which include both selection and distribution bias 
corrections. We find a linear relationship of the form 

A^2 = 0.0356(±0.0063) - 0.000012(±0.000001)c2:ls (22) 

where Afi 2 = Usft,corrected — Usfi IS the difference in distance 
moduli between the corrected and raw distances. The data 
and best fit line are displayed in Fig. 17. 

We can immediately see from (21) and (22) that there 
is a significant difference between the corrections. For small 
redshifts the Springob et al. (2007) correction is positive in¬ 
dicating raw distances are underestimated, while for large 
redshifts the correction is negative indicating raw distances 
are overestimated. Adjusting the intercept of Fig. 17(a) to 
zero we then find a hierarchy fj-sfi,corrected < Usfi < Ucf2 in 
the limit of large redshifts. This is consistent with the obser¬ 
vation of Watkins & Feldman (2015) that: “the distances are 
systematically larger in the Cosmicflows-2 catalogue [than 
in the COMPOSITE catalogue] due to a different approach 
to bias correction”. 

The fact that Fig. 22(b) has a negative slope means 
that the dominant correction cannot arise solely from se¬ 
lection and homogeneous distribution biases, since as noted 
above both of these effects underestimate true distances. 
The difference therefore must be due to the treatment of 
the inhomogeneous distribution bias, which Springob et al. 
(2007) have included but Tully et al. (2013) have not. 

The 1970 points in the SFI++ sample that are not also 
contained in the original CF2 survey, henceforth SFlH—hB, 
have been incorporated into CF2 without using their correc¬ 
tion (21). This data covers a larger range, up to redshifts of 
almost z = 0.1, whereas SFlH—HA only covers up to z = 0.06. 
Tully et al. (2013) state that these points, if corrected using 
(21) cause a “highly significant decrease in the Hubble pa¬ 
rameter with increasing velocity”. We independently verified 
this result^^. Thus Tully et al. (2013) do not adjust these 
distances, instead claiming that they are of a different nature 
altogether, the main difference being that these consist of 
cluster samples from a different survey (Dale et al. 1999a, b). 
In these samples rotation information for the galaxies was 
obtained from optical spectroscopy rather than the standard 
21cm Hydrogen line widths. However, Tully et al. state that 


so we are not interested in independently confirming this for our 
investigation. 

A decreasing Hubble constant below the scale of statistical ho¬ 
mogeneity is, to a limited extent, what is expected from the anal¬ 
ysis of Wiltshire et al. (2013). Thus trends which appear anoma¬ 
lous as compared to a standard FLRW expectation should not 
automatically be regarded as a signal of unaccounted observa¬ 
tional bias. However, there are also systematic differences that 
occur when binning in redshift, as in Figs. 18 and 19, as opposed 
to binning in distance with (4), so careful analysis is required to 
make sense of the different approaches. The direct calculation of 
the Hubble parameters in each bin is also different to that de¬ 
scribed in (4), as it is not clear which method Tully et al. (2013) 
use. We apply a simple weighted average of czi/r^ values and ob¬ 
tain consistent results (although as no values are tabulated by 
Tully et al. (2013) we can only verify by inspection). 
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Figure 17. The difference in distance modulus between: (a) the CF2 and raw SFI++ distances for the points common to the two 
surveys, as presented in Tully et al. (2013). (b) the SFIH—h raw and corrected catalogues for the same points. 


“it is not clear to us why this component of SFI++ does not 
manifest the selection Malmquist bias”. 

Since it appears two halves of the SFlH—h have been in¬ 
corporated into CF2 in different ways we determine whether 
they do actually show different characteristics. Considering 
A/i 2 for SFI-b-l-B we find a correction 

A/rz = 0.0417(±0.0061) - 0.000012(±0.000001)c2:Lg (23) 

which has an intercept consistent within uncertainties and 
identical slope to (22). Since there is no apparent difference 
using this test we repeat a similar analysis to that performed 
by Tully et al. (2013, Fig. 10) to test for bias. This analysis 
is based on the fact that selection bias is manifest by an 
increase in the Hubble parameter with redshift (Teerikorpi 
1993), for data binned by redshift. In Figs. 18 and 19 we 
repeat the analysis of the Hubble constant in redshift bins 
performed by Tully et al. (2013) for the subsets of interest, 
and compare the results. 

In Fig. 18 we produce plots equivalent to Tully et al. 
(2013, Fig. 10) for the SFI-|--|-A and SFI-|--|-B subsamples, 
both using raw distances. We subsequently find that the 
difference in the Hubble constant in individual redshift bins 
for the SFI-b-l-A and SFI-I—l-B ranges from 0.03 (j to 1.8 (t 
in individual bins. The weighted mean of these differences is 
0.84 (t, and thus we do not see a significant difference between 
SFH-bA and SFI-b+B. 

We note that two mutually consistent halves of the 
SFI-f-f sample have been incorporated into CF2 in different 
ways, but are uncertain as to what impact this may have on 
the full CF2 catalogue. 

On the other hand the COMPOSITE sample 
(Feldman et al. 2010) uses a much larger subset of SFI-I--I- 
distances corrected for Malmquist biases (after rejection of 
outliers). While it appears that there may be inconsisten¬ 
cies in the inclusion of the raw SFlH—h distances into the 
CF2 catalogue, it is possible that the SFI-f-l- corrected dis¬ 
tances are subject to systematic error also. It is for this 
reason that Springob et al. (2007) included both corrected 
and raw distances, to allow others to take on the challenging 
task of Malmquist bias corrections. However, the CF2 cat¬ 
alogue is only corrected for selection bias, leaving the cor¬ 


rection of the distribution homogeneous and inhomogeneous 
Malmquist biases as a task for the user. 

As another test of differences between the SFI-f-fA and 
SFI-b-f-B subsamples we have also repeated the analysis of 
Fig. 18, but now to compare the raw and corrected distances 
within each subsample. Fig. 19 shows the comparison for 
the SFIH—hB subsample. In this case the difference in the 
Hubble parameters in each bin vary from a minimum O.Olcr 
to a maximum 2.2a in individual bins, with a weighted mean 
difference of la. For SFI-f-fA the difference in the Hubble 
parameters vary from O.OSa to 1.9a in individual bins, with 
a weighted mean difference of l.Ocr. Thus again we do not 
see a significant difference between the subsamples. 

It may appear surprising that the raw and uncorrected 
data only differ by Icr on average when binned by redshift. 
However, once inhomogeneous Malmquist bias is accounted 
for the sign of the correction is different at large redshifts as 
compared to low redshifts, meaning that for an intermediate 
range the correction is small. The approach by Tully et al. 
(2013) of binning in redshift is not an appropriate one to use 
when performing a parameter minimization that involves 
boosts to rest frames in which the redshift is changed, as in 
§3-4. Rather we followed Wiltshire et al. (2013) in binning 
by distance. This led to differences between the raw and cor¬ 
rected data sets which are statistically much more marked 
than is evident if one bins by redshift. 


7 DISCUSSION 

We have investigated the extent to which it is possible to 
define a standard of cosmological rest based on a frame in 
which variation of the spherically averaged (monopole) non¬ 
linear Hubble expansion is a minimum. Such averages do not 
make any assumptions about the geometry of space below 
the scale of statistical homogeneity (< 100h“^Mpc) of the 
sort implicit in the standard peculiar velocities framework, 
which assumes a Euclidean spatial geometry. 

We studied the systematic variation that arises when an 
arbitrary boost is made from a frame of reference in which 
the spherically averaged Hubble expansion is most uniform. 
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Figure 18. The Hubble parameter, Hi = czijri, computed for each individual data point (coloured points) and from averaging in 
1000 km sec“^ bins (black points) using the: (a) SFIH—|-A subsample; (b) SFIH—hB subsample, (c) Comparison of the averaged points 
in (a) and (b) with blue crosses being from SFIH—hA and black filled circles from SFIH—hB. 
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Figure 19. The Hubble parameter, Hi = czi/ri, computed from 
averaging in 1000km sec“^ bins using the SFIH—hB subsample 
with corrections for Malmquist bias (blue crosses) and without 
corrections (black filled circles). 


We found that such a systematic variation is indeed de¬ 
tected to a statistically significant extent between the CMB 
and LG frames of reference, using the COMPOSITE sample 
(Watkins et al. 2009; Feldman et al. 2010). 

This supports the proposition made by Wiltshire et al. 
(2013) that the Local Group frame may be a more suitable 
cosmic rest frame than the standard CMB frame, and that 
consequently a significant fraction of the CMB dipole may be 
nonkinematic. This proposition can be directly tested in nu¬ 
merical ray-tracing simulations (Bolejko et al. 2015). Given 
that the CMB frame is still the de facto choice for the cosmic 
rest frame, this conclusion would have a far reaching impact 
for many areas of cosmology (Wiltshire et al. 2013), includ¬ 
ing the question of large angle CMB anomalies (Ade et al. 
2014b). 

We extended the analysis to search for an improvement 
on the Local Group as the standard of rest, using a vari¬ 
ety of tests. Very large boosts from the LG frame can be 
excluded if we simultaneously demand that while the resid¬ 
ual variation of the Hubble expansion should be small in 
the regime of nonlinear expansion at small distances, a clear 
signature of an emerging linear Hubble law should also be 
found at larger distances. However, we found that in apply¬ 
ing all possible tests there is still freedom to perform quite 
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large boosts close to the plane of the galaxy, presumably 
because the lack of data in the Zone of Avoidance leads to a 
lack of suitable constraints there. This hypothesis could po¬ 
tentially be tested on simulated data using exact solutions 
of Einstein’s equations (Bolejko et al. 2015). 

Since our conclusions depend on the COMPOSITE 
sample, we repeated our analysis insofar as it was possible 
for the recently released Cosmic flow s-2 sample (Tully et al. 
2013). This catalogue of 8162 galaxy, groups and cluster dis¬ 
tances is considerably larger than the COMPOSITE sample, 
and potentially could provide more accurate results although 
it is also of course limited in the Zone of Avoidance. 

We have found very significant differences in the re¬ 
sults for the CF2 and COMPOSITE samples, as can be 
seen by comparing Fig. 3 of Wiltshire et al. (2013) and our 
Fig. 11. These result from differences in the treatment of the 
Malmquist bias between the SFI++ catalogue and CF2 cat¬ 
alogue, as previously noted by Watkins & Feldman (2015). 
We also found apparent inconsistencies in the manner of 
inclusion of subsamples of the SFI-f-|- catalogue into the 
CF2 catalogue, with respect to the treatment of Malmquist 
bias. More significantly, since the reported CF2 distances 
do not include corrections for the inhomogeneous distribu¬ 
tion Malmquist bias they may be of limited use until such 
corrections are applied. 

The conclusions of Wiltshire et al. (2013) are dependent 
on the treatment of the Malmquist bias in the SFI-I—I- cat¬ 
alogue being accurate, as this constitutes the largest part 
of the COMPOSITE sample. Naturally, one might ques¬ 
tion whether any systematic procedure of Springob et al. 
(2007) could somehow spuriously lead to an unusually uni¬ 
form Hubble expansion in the LG frame through an error in 
the Malmquist bias procedure. 

We find no grounds for this. In particular, our analysis 
shows that the difference between the CMB and LG frames 
has the distinctive signature of a systematic boost offset (12) 
noted by Wiltshire et al. (2013). Nothing in the Malmquist 
bias correction procedure of Springob et al. (2007) could ob¬ 
viously introduce this signature through a systematic error. 
Their analysis does not single out the LG, or LS, frame in 
any way; indeed all their redshifts are referred to the CMB 
frame. Furthermore, although the remaining bias means that 
the CF2 sample is currently unusable from the point of view 
of determining a frame of minimum spherically averaged 
Hubble expansion variation in the nonlinear regime, we ob¬ 
served in §5.2 that the difference of the CMB and LG frame 
spherical averages nonetheless still shows the signature of 
the systematic boost offset in the CF2 data. 

Since the boost offset is detectable in the independently 
reduced CF2 data, it cannot be an artefact of the Malmquist 
bias treatment of Springob et al. (2007). Furthermore, the 
departure of the nonlinear expansion from the boost power 
law (15) that is seen when comparing Figs. 2 and 14 is pre¬ 
cisely what is to be expected if there are additional unac¬ 
counted uncertainties in individual CF2 distances as com¬ 
pared to the SFI-I—I- ones: the distance range of structures 
associated with the nonlinear expansion is broadened. 

If Fig. 11 was based on accurate distances, it would im¬ 
ply that the Hubble expansion in all frames of reference is 
far less uniform than might reasonably be expected in any 
viable cosmological model; in particular, there is a monopole 
or “Hubble bubble” variation of order 15-20% in the range 


20 < r < 60 h“^Mpc in the CF2 sample, as compared to 4- 
5% in the GOMPOSITE sample. The largest “Hubble bub¬ 
ble” variation that has ever been claimed on such scales 
using more accurate Type la supernovae distances in the 
CMB frame is 6.5 ± 2.2% (Zehavi et al. 1998). 

Tully et al. (2013) chose not to correct for the dis¬ 
tribution biases, as they wished to separate “the issues 
of distance measurements and velocity field inferences”. 
Indeed, in the peculiar velocities approach the distribu¬ 
tion bias may be much less significant. In new work 
Hoffman, Courtois & Tully (2015) use the CF2 catalogue to 
reconstruct large scale structure by means of the Wiener fil¬ 
ter and constrained realizations of Gaussian fields assuming 
a WMAP constrained AGDM model as a Bayesian prior. 
They observe that “the Malmquist bias introduces a spu¬ 
rious strong monopole term into the reconstructed veloc¬ 
ity field but is expected to hardly affect the bulk velocity 
which is associated with the dipole of the velocity field” 
(Hoffman et al. 2015). This would appear to be the coun¬ 
terpart of the large monopole we observe in Fig. 11 in our 
analysis. Hoffman et al. (2015) corrected for the bias but 
noted that the bulk velocity analysis is “virtually unaffected 
by the Malmquist bias”. 

In our case, residual bulk flows on scales > 100 /i“^Mpc 
may be an artefact of using the GMB rest frame as the stan¬ 
dard when it does not coincide with the frame of minimum 
Hubble expansion variation. In particular, Wiltshire et al. 
(2013) observed that a dipole fit to the Hubble expansion 
gave a dipole magnitude which increased to a residual value 
in the outer shells as the shell radius was increased in the 
CMB frame, but not in the LG frame. Thus large scale bulk 
flows are not our primary interest. Rather, we are interested 
in detecting the systematic monopole variation (12). In dis¬ 
tinction to the peculiar velocity approach our method by ne¬ 
cessity is sensitive to a monopole bias. In fact, our method 
of binning in radial shells by distance with an anchoring to 
Hq, is particularly sensitive to any distribution bias which 
follows from a number density, N{r), with strong gradients. 
The bias effect in Fig. 11 can be largely reproduced by apply¬ 
ing a uniform Hubble law to the CF2 redshifts, adding Gaus¬ 
sian scatter to create a mock distance catalogue, and then 
applying our binning strategy (R. Watkins, private commu¬ 
nication) . 

It may be possible to construct the 635 km sec“^ ve¬ 
locity attributed to the LG within the AGDM model, as 
has recently been claimed by Hess & Kitaura (2016) who 
used constrained A-body simulations and nonlinear phase 
space reconstructions to arrive at a value ~ ^ 

137km sec“^. However, this itself does not constitute a 
proof of the standard kinematic interpretation, but rather a 
verihcation within the 20% uncertainty of a computer sim¬ 
ulation. The ACDM model is certainly phenomenologically 
very successful, and any competing model can only be viable 
insofar as many of its predictions are close to the standard 
model, as is the case, for example, in the timescape cosmol¬ 
ogy (Wiltshire 2009; Smale & Wiltshire 2011; Smale 2011; 
Nazer & Wiltshire 2015). What is important in testing the 
standard model is to seek observations which are not ex¬ 
pected in its framework. Although the signature of the sys¬ 
tematic monopole boost offset (12) between the GMB and 
LG frames should be checked in ACDM simulations, it is not 
an observation that should obviously arise if all cosmologi- 
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cal redshifts arise purely from a FLRW geometry plus local 
Lorentz boosts at any scale. Properly characterizing and de¬ 
termining a frame of minimum nonlinear Hubble expansion 
variation is therefore a fundamental question open to more 
precise observational tests in future. 

In conclusion, if we make no assumptions about our own 
standard of rest, but make arbitrary boosts at our location, 
then analysis of the COMPOSITE sample yields a degener¬ 
ate set of frames for which the variation of the spherically 
averaged nonlinear Hubble expansion is minimized. This set 
of candidate cosmic rest frames includes the LG frame but 
excludes the CMB frame, as was already established with 
very strong Bayesian evidence by Wiltshire et al. (2013). 
The degeneracy in defining the minimum variation frame 
is associated with a freedom to perform boosts in the plane 
of the galaxy without changing the statistical significance of 
the fit, probably due to the lack of constraining data in the 
Zone of Avoidance. 

The larger CF2 sample may potentially tighten the con¬ 
straints on the dehnition of the frame of minimum nonlinear 
expansion variation. However, it is first necessary to reduce 
the data in the manner of Springob et al. (2007) to remove 
the inhomogeneous distribution bias which appears to be 
the source of the large discrepancies found in §5. The Uni¬ 
verse is very inhomogeneous below the scale of statistical 
homogeneity, and ironically it is only once the biases associ¬ 
ated with such inhomogeneities are removed that a picture 
of a remarkably uniform average Hubble expansion actually 
emerges. 

A careful treatment of inhomogeneous Malmquist bias 
is therefore key for the future progress of our understanding 
of the nature of cosmic expansion as the surveys grow ever 
larger. 
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APPENDIX A: LINEAR REGRESSION VIA 
TOTAL LEAST SQUARES 

Consider a general linear model with errors in both the de¬ 
pendent and independent variables. We can express such a 
model as 

yi = Po + hixi 

{Yi,Xi) = {yi,Xi) + {ei,Ui) 

where (U, Xi) are the observations, {yi,Xi) are the true val¬ 
ues and (ei,Ui) are the measurement errors. We assume the 
measurement errors to be normally distributed with a co- 
variance matrix 



We can extend this analysis to allow for different errors 
at each point, which we will refer to as weights given by 
u}{Xi) = l/auu,i, oj{Yi) = and the correlation coeffi¬ 

cient between the errors given by 7 i = aeu,i 

To carry out a least squares minimization we must cal¬ 
culate the statistical distance from an observation to the true 
value. In standard least squares this distance is the vertical 
distance from the data point to the model, since the values 
of the independent variable are assumed to be exact. Now, in 
the most simple case we would have the squared Euclidean 
distance from the observed data points to true value in the 
model as 

[Yi - (/3o -b Pixi)^ + {Xi - Xif = 4 + ul (A3) 

but if the variances of d and Ui are different from unity 
this statistical distance becomes el + ^ and if these 

variances are correlated we must use the covariance matrix 
to give the “statistical” distance 

{Yi - 00 - PiXi, Xi - Xi) {Yi -00- 0iXi, Xi - xW ■ 

(A4) 

In our case we must determine the values of 0o and 0i that 
minimize (A4). That is, we must find the values {xi, yi) and 
{ 00 , 0i) that minimize this sum for the given observations. 


First we fix the Xi values by treating them as unknown con¬ 
stants in a standard linear regression of the form 


Yi-00 

Xi 


01 

1 


Xi -b 


Ci 

Ui 


(AS) 


for which the generalized least squares estimator gives 
Xi = [{0i,l)Y-40i,W] {0i,l)Y-4Yi - 0o,XW ■ (A6) 

Substitution of Xi into (A4) gives 

(W -00- 0,X,.f 

(rTee ‘^011Xeu “b 0\(Xuu) 

SO that after summing over all N points we obtain 

o Eti -ho-piXi)^ ..Si 

(o-ee — 2/3lCTeu + 0i(Xuu) ’ 

which is the expression to be minimized. In (York 1969) the 
linear equation that minimizes (A8) is given by 


01 = 


where 


Ef=i 

r Ui , dVi 7iYi 

[uiiYi) ' ui(Xi) Ci \ 



Ui 1 PiVi PiliUA 

iij(Xi) oci J 


(A9) 


Ki =uj{Xi)uj{Yi), Ui = Xi-X, Vi = Yi-Y, 

/ N \ ^ / N \ N 

X = iY^ZiXi] /Y,Xi and Y = ZiYi] 


Zi = 


w{Xi)w{Yi) 


uj{Xi) -b b^uj{Yi) — 2b^iai' 


Clearly (A9) requires an iterative process to find 0i which 
begins with an initial guess which may be found from per¬ 
forming a standard linear regression assuming the Xi to be 
exact. After 0i is obtained the value of 0o is found from the 
fact that the mean must be on the best fit line and thus 
00 = Y- 0iX. 

The uncertainties in the parameter values, and cr/jj, 
are (Titterington & Halliday 1979) 


_ T,ZWi _ 

ij:zix^i){j:zi)-{j:ziXi)^’ 

YI 

ij:zix^i){j:zi)-{j:ziXi)^- 


(AlO) 

(All) 


We now return to the transformed model from (12) 
which takes on the form 


log{SHi) =plog((r?)) -b log 



(A12) 


such that we may identify yi = log{SHi), Xi = log((rf)), 
and (^i,^o) = (^p,log (^ 2 ^^^=(p,logA). 

Since our distances are given in units /i“^Mpc, indepen¬ 
dent of the overall normalization of Hq, rather than directly 
performing a fit to (15), (A12) we instead perform a fit to 
the relation 


{H{ - Hs) 



(A13) 


This gives a direct estimate of v in km sec from which 
we obtain A ~ /{2Hg)- 
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ERRATUM 


The paper “Defining the frame of minimum nonlinear 
Hubble expansion variation” was originally published in 
MNRAS, 457, 3285 (2016). 

There is a small numerical error^® in the denominator 
of the final line of eq. (12) of McKay & Wiltshire (2016) 
where the factor “2” should be “3”, so that it reads^® 


AH, = H',-Hs 


E 

i=l 


{v cos (fliY 


i^i)s 



( 12 ) 


Likewise three lines below eq. (13), we should have 
{{v cos cj)ifi) ^ Jf sin while in the 

sentences below eq. (15) and eq. (A13) we should have 
A « ii^/(377j,), and eq. (A13) itself should read 


Ho {H's ((r?)^ 


(A13) 


Correction of this small error does not affect any of 
the conclusions of the paper, but it does increase the value 
of the boost velocity, v, inferred from the parameter A in 
eq. (12) or (A13) by 22.5%. This has no material conse¬ 
quence, however, since the uncertainties in v are very large. 
In particular, using the best-fitting value of A, in place of 
the values given in section 3 we find v = Qdfilggg® km sec“^ 
(v = 5231 km sec“^) for the primed (unprimed) shells re¬ 
spectively, both still being consistent with the actual boost 
magnitude of 635 ± 38 km sec“^. 

The only instance in the paper where incorrect derived 
values of the boost velocity were used in any further nu¬ 
merical analysis was in Section 4.3, where these values were 
assigned to Uder and compared to the true boost velocity 
Wtrue, as parameters were varied in the {v,£,b} parameter 
space. The 22.5% correction to Uder results in slight changes 
to Fig. 9(a), as shown, but no significant change to Fig. 9(b). 

Using the corrected value to define the frame X of sec¬ 
tion 4.3, we find that Uder/Utrue = 1 now occurs at Vtme = 
239km sec“^ in a direction (£,b) = (59°,— 3°) ± (5°,4°), 
which is only changed by 1°. For this boost p = —1.25 ± 
(0.92)stat ± (0.37)sys, which is still consistent with p = —1 
within the large uncertainties. (With only the primed shells 
we now obtain p = —1.17 ± (1.06)stat.) 

None of the physical conclusions concerning frame X 
are changed from Section 4.3 of McKay & Wiltshire (2016). 
In particular, the changes to Fig. 10(b),(c) are so small 
that they are indiscernible apart from a very small offset 
of the first data point. Moreover, as before frame X is only 
weakly disfavoured compared to the global Xa minimum, 
ln(Pj^.y/P^) = 1.7, and is within 1 cr of the global minimum 
of Xbi similar to the LG frame on both counts. 



Figure 9(a) Corrected. Ratio of the derived velocity from the 
best fit power law to the true boost velocity. At each v the direc¬ 
tion is determined by the best fit for which p is closest to —1. 
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The error was first made in eq. (9) of Wiltshire et al. (2013), 
and was later also repeated in eq. (2.8) of Bolejko et al. (2016). 
However, as the value of v inferred is not used in either of these 
papers, none of their results are affected in any way. 

This relation has now been extended by Kraljic & Sarkar 
(2016) to include the case that a bulk flow exists in the frame 
of minimum spherically averaged Hubble expansion variation. 
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